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Abstract 

We study the two-dimensional lattice Gross-Neveu model with Wilson 
twisted mass fermions in order to explore the phase structure in this setup. In 
particular, we investigate the behaviour of the phase transitions found earlier 
with standard Wilson fermions as a function of the twisted mass parameter 
/X. We find that qualitatively the dependence of the phase transitions on is 
very similar to the case of lattice QCD. 
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1 Introduction 



Recent results from simulations of lattice QCD with Wilson twisted mass fermions 
m 12] revealed an intriguing phase structure with an Aoki phase at strong coupling 
jSl E] and pronounced first order phase transitions at intermediate values of the 
coupling which weakens when the continuum limit is approached jSJ IHl El E] • These 
findings are supported by calculations in chiral perturbation theory for standard 
Wilson fermions [H] which were later extended to the case of Wilson twisted mass 
fermions [TUl CH [121 UHl UH UH] • See refs. [121 [I7[ for reviews. The appearance of the 
first order phase transition turns out to be a generic phenomenon of Wilson fermions 
and takes place independent of whether the twisted mass parameter is switched on 
or not. 

Of course, one may argue that on a finite lattice no phase transition can occur 
and that the metastabilities observed in refs. [SI [El [3 [HI are purely an artefact of the 
algorithm employed. From the analytical side, chiral perturbation theory (xPT) 
is a powerful tool to determine the phase structure. However, the assumption of 
xPT is that aA ^ 1, where a is the lattice spacing and A is the scale parameter. 
Moreover, the analysis by xPT is done up to 0{a?) only. Thus, it is unclear whether 
at large values of a and by inclusion of higher order terms the picture can change 
substantially. 

Motivated by these considerations and in order to have another angle on the 
problem, we decided to study the phase structure of the Gross-Neveu model [TH] 
on the lattice including a twisted mass term, taking this as a toy model of twisted 
mass lattice QCD. In the large-A^ limit, where is the number of flavours, the 
Gross-Neveu (GN) model is asymptotically free and has bound states (TH]. The 
GN-model is one-loop exact in the large-A^ limit and hence, in this limit, higher 
order contributions of the lattice spacing are taken into account. In this sense, the 
large- A^ investigation of the GN-model goes beyond xPT and can provide therefore 
a complementary picture of the phase diagram, although one should always keep in 
mind that we are considering a two-dimensional model only. 

The GN-model has been studied for pure Wilson fermions in the large- A^ limit 
in refs. |19| [^. In ref. [2I|, the phase structure of the lattice GN-model with Wil- 
son fermions is first predicted and a so-called Parity- Flavour broken phase (or Aoki 
phase) was found. A refined analysis fHl showed that even the zero temperature 
phase diagram is rather complex and reveals besides the second order phase tran- 
sition to the Aoki phase strong first order phase transitions. In this work we want 
to extend the investigations of refs. [2B [201 ^| to the case when a twisted mass 
parameter /i is switched on. In particular, we want to address the question of the 
/^-dependence of the first order phase transition. We hope that our findings can give 
some, at least qualitative insight also into the generic phase diagram of lattice QCD 
with Wilson fermions. 
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2 Lattice Gross— Neveu model with twisted mass 
fermion 



A formulation of the two-dimensional Gross-Neveu model on the lattice with Wilson 
fermions is given in refs. [211 1201 ^] • Our study here follows refs. fIU\ IT^ where the 
the lattice couphngs of the model are treated separately. We choose the so-called 
twisted basis [21 for our investigation in which the Lagrangian of the model reads 

2 2 

(1) 

where (x, x) denote the fermion fields in the "twisted basis" or x-basis, and is the 
number of the "flavours" . The parameters of the model are the critical fermion mass 
Mcr, the off-set fermion mass m^, the twist angle oj and the couplings g„ and Qt,. In 
this version of twisted mass fermions, the Wilson term is not twisted. An important 
remark is that in general the coupling constant 7^ Q-k on the lattice PUIITU] . 
In momentum space, the Lagrangian in the x-basis is given as 

2 2 

Cgn = X [K{p) + W{p) + mo + 275r3/i] X - ^{xxf - ^(x^TsTsx)', (2) 

where K{p) = z X],/ 7^ sin(p,,a), W{p) = rJ2ui^ ~ cos(p^a)), mo = Mcr + mgcoscj, 
/i = rriq sin uj and r the Wilson parameter. The situation of "full twist" is realized 
when u = ±7r/2 or, equivalently, mg = M^r- 

In the large- A^ limit we obtain the Lagrangian 

Cgn = X [Kip) + W{p) + ctl + 275X37?^] X + ^{(^L - mof + ^{ttl - (3) 

where, by the equations of motion, the bosonic fields are defined as 

(^L = ~^{XX) , 711 = -^(x^75i"3X) + /i • (4) 

To solve the model, we integrate out the fermion fields, treating the vr- and a- 
fields as constant since the quantum fluctuations are suppressed in the large- A^ limit. 
We then obtain the effective potential Ve//, 

Z = J [dx] [dx] exp{- J d\CGN} = e^p{-NnVeff} , (5) 
where Q is the space-time volume. 

Veff in large- A^ limit is given by (setting here and in the following r = 1), 

Veff = ^{(TL-mof + ^{til- fif (6) 
- / ln[^ sin^ip^a) + {^(1 - cos(p^a))}2 + {al + nl) + 2(Tl ^(1 - cos(p^a))], 

•J'O 
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where = J^^ "(l^p^- Note that eq. (jHl) has the same form as the effective potential 
in the Wilson fermion formulation of the GN-model on the lattice |2()| ITHj . The 
difference shows up only in the definition of the field ttl, where there is an explicit 
appearance of the twisted mass parameter fi. 



2.1 Phase structure and Pion mass in twisted basis 

At /i = 0, i.e. for pure Wilson fermions, the phase structure of the Gross-Neveu 
model is already quite complicated showing first order and second order phase tran- 
sitions. It is unclear, whether the presence of the twisted mass term can change the 
phase structure. In order to explore the phase structure of the lattice GN-model, 
it is sufficient to solve the saddle point equations as they can be derived from the 
effective potential Ve//. 

In the following we will concentrate on the twisted basis, because then the effec- 
tive potential has the same form as the one for ordinary Wilson fermion. The effect 
of the twisted mass parameter is, to 1-loop order, only included via the field hl by 
a shift. 

The saddle point equations are 

-(rLF{aL,nL)-G{(rL,nL)=0, (7) 

TTlFK, TTi) = 0, (8) 



dVeff 


aL -niQ 


dcTL 


9l 


dVeff 




d-KL 


9l 




r 2 


J 


Ip A{aL,7rL) 


= ( 





where 



2Ei.(l - cos(p^a)) 



-2(Tl^(1 - cos(pi.a)) . 



In the saddle point equations the 1-loop tadpole contributions are canceled by 
the potential terms of and tt^. The solution of the saddle point equations is 
non-trivial and will differ from the one obtained with pure Wilson fermions opening 
thus the possibility of new phenomena. 

The neutral pion mass and the sigma particle mass are computed from the second 
derivative of the effective potential at the solution of the saddle point equations. 
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Our strategy is the following. We first fix the values the couplings g1 and g1 
such that the perturbative tuning relation between g"^ and gl. as computed for the 
ordinary Wilson fermion formulation holds, 

^ = 4 + 4C2 + 0(a), C2 = ^ + -^. (11) 
gi g^ 27 127r 

In the following table, the (perturbative) tuning values for various g^ are given. 



9i 


0.1 


0.5 


1.0 


2.0 


9l 


0.094 


0.381 


0.617 


0.893 



We then also fix the value of the twisted mass /i and determine the values of 
and TiL from the solutions of the saddle point equations eqs. ((71 |HI) as a function of 
rriQ. In this way, we can study whether ctl or til undergo a phase transition. Once 
we have the value for ctl and vr^ we can insert them into the gap equations, eqs. (0 
fTI]|) to compute also the corresponding masses. 

Let us remark that in principle there is also an additional tuning relation between 
the coupling gl and the fermion mass mo. To explore the phase structure, we will, 
however, vary mo and hence will not make use of this tuning relation. We will also 
restrict ourselves to values oi g^ < 2, since for larger values of the coupling the 
tuning relation of eq. ^\\\ will not hold anymore. 



3 Results 

In this section we will give our results obtained from the solution of the saddle point 
equations. We used a lattice of size 128^ and solved the equations numerically. We 
always checked that we indeed found the correct solutions by inspecting the extrema 
of the effective potential itself. 

In fig.Qwe show the expectation value (ttl) of the vr-field as a function of the mass 
parameter mo for various values of the twisted mass parameter /i. Let us summarize 
the results of our large- computations. At /i = and increasing the values of 
mo, we find in accordance with ref. first and second order phase transitions, 
where the vr-field assumes a non-zero value. For larger values of mo there is a first 
order transition where (vtl) jumps back to a zero value. This structure changes in 
two aspects when the twisted mass parameter /i is switched on. First, the second 
order phase transition vanishes and (vtl) shows an analytical behaviour. Second, 
for some range of < /icrit the first order phase transition remains. However, it 
becomes weaker when /i is increased and turns into a second order phase transition 
when a critical value /icrit is reached. For /i > /icrit the phase transition disappears 
altogether and turns into a completely smooth and analytical behaviour. 
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When the couphng g1 is made smaller, with the associated tuning of the 
region in mo where the phase transitions occur shrinks and also the jumps at the 
first order phase transitions become smaller. Thus the qualitative picture does not 
change when we approach the continuum limit. Note that also the value of /icrit gets 
smaller with decreasing coupling g^. 

The fate of the first order phase transition is very much reminiscent of the 
situation of four-dimensional lattice QCD. Also there it has been found from nu- 
merical simulations [3 El E] and from calculations in chiral perturbation theory 
jni El Oni El E| that the first order phase transition terminates at a critical values 
of /i. Note that the values of mo where the first order phase transition takes place 
moves to larger values when /i is increased. A similar phenomenon can be expected 
also in lattice QCD. What is very interesting from our analysis in the GN-model 
-and apparently different from current findings of lattice QCD and 0(a^)-xPT- 
is that this first order phase transition is accompanied by a second order phase 
transition at = or its remnant at // > 0. 

In fig. 121 we show at the example of the coupling values gl. = 2, g^ = 0.893 the 
behaviour of the pion mass (top), the a-field expectation value (middle) and the 
(T-mass (bottom). At /i = and increasing mo, the pion mass first goes to zero 
and vanishes at the second order phase transition. It then assumes non-zero values 
again until it jumps when the first order phase transition is hit. Switching on fi, the 
second order phase transition vanishes and the behaviour of the pion mass becomes 
continuous. The first order phase transition remains, however, until fi = /icrit- For 
larger values of n also the first order phase transition vanishes and the pion mass 
shows an overall analytical behaviour. 

The cT-field expectation value and mass do not show a large sensitivity on the 
second order phase transition. However, this field feels the first order phase transi- 
tion strongly and shows jumps in the field expectation value as well as in the mass. 
The picture shown here for the couplings gl = 2, gl = 0.893 stays qualitatively 
unchanged when the coupling gl is reduced with the associated tuning of gl- 

Finally and for completeness, we show in fig. El the fate of the so-called a-phase 
transition ^U] for a value of the coupling gl = 2. Note that the value of the coupling 
is now gl = 0.8 and does not correspond to the ones of the perturbative tuning. This 
phase transition is characterized by the fact that at /i = the cr-field expectation 
value ((Tl) shows a discontinuous jump, whereas the vr-field expectation value stays 
essentially zero. This phase transition survives when /i is switched on and can clearly 
be seen in the jump of ((T^). It is interesting to observe that the point in mo where 
this jumps happens is moved to larger values of mo when fi is increased. However, 
also for this case the phase transition vanishes when /i is increased above a certain 
critical value. Furthermore, the vr-field shows now at /i > also a jump and exhibits 
the first order phase transition. 
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4 Summary 



Following the earlier work of refs. j2T| l2(H IT9] we studied in this letter, how the 
presence of the twisted mass term influences the known phase diagram of the two- 
dimensional Gross-Neveu model which we consider as a toy model of lattice QCD. 
We found that the coexistence of the Aoki phase and the first order phase transition 
as it is obtained at vanishing twisted mass parameter /x ^H] is not kept at /i > 0. The 
first change is that at /i > the second order phase transition turns into a smooth 
and analytical behaviour, although e.g. the vr-field expectation (vr^) still shows a 
rapid change. The second change is that the gap in (vr^) of the first order phase 
transition shrinks with increasing values of fi and finally vanishes at a critical value 
of = Atcrit- In addition, the values of mo, where the first order phase transition 
takes place, is shifted to larger values when fi is increased. 

The motivation of our investigation has been to provide a complementary anal- 
ysis of a "QCD-like" model to the findings of numerical simulations El 13 E] and 
from (Wilson) chiral perturbation theory P [HIl EH 1121 CHI 113 EH! • Both of these 
approaches have some shortcomings: in the numerical simulations it is not clear 
whether the metastability effects are due to a bad behaviour of the algorithms em- 
ployed. xPT on the other hand can, in principle, only predict the phase structure 
reliably in a region close to the continuum limit where the O(a^) expansion is valid 
and the quark masses are small. 

Our findings from the large- iV analysis of the GN-model as described above show 
a striking similarity to the results from the four-dimensional simulations and from 
xPT. This includes the shift in the critical mass, the weakening of the first order 
phase transition as function of increasing /i and the vanishing of the first order 
phase transition altogether when /i > /icrit- However, there is also a very noticeable 
and important difference to the picture we have presently in four-dimensions: this 
concerns the coexistence of the first order phase transition with a second order 
transition to the Aoki phase. It is unclear to us whether this difference is an artefact 
of the GN-model in the large- analysis. Another, very interesting possibility is of 
course that so far both, the numerical simulations and the analysis in xPT have 
missed this possibility. Clearly, if the Aoki phase transition is very close to the first 
order phase transition it will be very difficult for the numerical simulations to resolve 
both. It could be speculated that the coexistence of both phases will show up in 
xPT only when higher orders are taken into account*. From our analysis here we 
would just like to give a warning that maybe the phase structure of Wilson lattice 
QCD is even more complicated than our present picture suggests. If the picture 
of the phase diagram from our analysis in the GN-model would be correct also for 
lattice QCD, it would provide an argument for using a no n- vanishing twisted mass 
parameter, since then one would avoid the complication of the existence of two phase 

*Indeed, in some recent calculations in xPT for four-dimcnsional latticc-QCD, taking higher 
order effects into account such a coexistence of the Aoki phase and the first order phase transition 
-as we found here- seems to take place |22| . 
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transitions. 

We also want to mention that the a-phase transition, which is present in the 
GN-model with Wilson fermion, changes to a normal first order phase transition, 
when n is switched on. 

In this letter we did not address a number of additional interesting points that 
would deserve a further investigation: in order to determine the phase transition, in 
particular in order to study the second order phase transition, the second moment, 
for instance the susceptibility, may be better. In principle one could also study the 
scaling behaviour of physical quantities analytically toward the continuum limit. In 
addition, an analysis within perturbation theory could provide important insight 
into the renormalization properties of the GN-model with the twisted mass term. 

Finally, it would be straightforward to change the discretization of the fermions 
in the GN-model. For example, one could study a chiral invariant GN-model. In 
this way, it would be possible to learn how the change of the fermion part of the 
action changes the phase structure. This is clearly important for four- dimensional 
lattice QCD simulations. 
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Figure 1: (tt/,) as a function of the fermion mass mo- The couphng g^^ is the per- 
turbatively tuned one obtained from eq. (fTT|l and assumes values of gl = 0.893, 
gl = 0.617 and gl = 0.381 from top to bottom. We represent some values of the 
twisted mass parameter fi by open symbols. Other, intermediate values of /i are 
represented by solid lines in the graph. 
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Figure 3: The vr- (top) and a-field (bottom) expectation values as a function of mo 
at the so-called cx-phase transition [T^ . 



11 



